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Broken symmetry states characterizing density waves of higher angular momentum in 
correlated electronic systems are intriguing objects. In the scheme of characterization by 
angular momentum, conventional charge and spin density waves correspond to zero angu¬ 
lar momentum. Here we explore a class of exotic density wave states that have topological 
properties observed in recently discovered topological insulators. These rich topological 
density wave states deserve closer attention in not only high temperature superconduc¬ 
tors but in other correlated electron states, as in heavy fermions, of which an explicit 
example will be discussed. The state discussed has non-trivial charge 2e skyrmionic spin 
texture. These skyrmions can condense into a charged superfluid. Alternately, they can 
fractionalize into merons and anti-merons. The fractionalized particles that are confined 
in skyrmions in the insulating phase, can emerge at a deconflned quantum critical point, 
which separates the insulating and the superconducting phases. These fractional par¬ 
ticles form a two-component spin-singlet chiral {d^ 2 _y 2 di idxy) wave superconducting 
state that breaks time reversal symmetry. Possible connections of this exotic order to the 
superconducting state in the heavy-fermion material URu 2 Si 2 are suggested. The direct 
evidence of such a chiral superconducting state is polar Kerr effect that was observed 
recently. 

Keywords: Skyrmion; Merons and antimerons; Density-wave; Chiral superconductivity; 
Deconflned quantum criticality. 


1. Introduction 

In this review we address novel role that skyrmionJ^ and their fractionalized avatars, 
merons and anti-merons, can play in certain correlated electron systems. This re¬ 
view focuses on our own worlP® except where contributions by other authors have 
provided ingredients. The exploration of skyrmions in the density wave states of 
the kind discussed here has, to our knowledge, not been discussed in places other 
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than in our own work. One of the key features is that it can predict a chiral d-wave 
superconducting state that breaks time reversal symmetry (TRS), which appears to 
have been directly observed in recent polar Kerr effect (PKE) measurements by the 
Stanford groupPin a heavy fermion material: URu 2 Si 2 . Broken TRS is a necessary 
condition for a nonzero PKE.I^ 

The order parameter in URu 2 Si 2 (URS), a heavy-fermion material, below the so 
called hidden-order (HO) transition at Tho = 17-5 K is unknown despite its discov¬ 
ery over a quarter century ago; see Ref. [Sand Ref. |8|and references therein. Within 
this phase lies a much less explored unconventional superconducting state with a 
transition temperature Tc 1-5 kEIHIII It stands to reason that there must be an 
intimate relation between the two. While numerous theoretical models have been 
proposed to explain the HO phasepHU there are very few attempts to explain the 
mechanism of the unconventional superconductivity. It is our central interest to ex¬ 
plore the connection between the two states to provide a skyrmionic mechanism for 
the unconventional superconducting state, which arises from an intriguing density 
wave state, termed mixed singlet-triplet d-density wave (st-DDW).l2lElThis state has 
no net charge or spin modulations and does not break TRS. It does have topological 
order with quantized spin Hall effect.^ Thus, it is naturally impervious to common 
experimental probes and can be aptly described as an hidden order state. Determi¬ 
nation of the density wave state posited here may be possible through two-magnon 
Raman scattering, nuclear quadrupolar resonance, or the skyrmions themselves. In 
a more general context, our work reflects the rich possibilities of emergent behavior 
in condensed matter systems. Attempts were made to describe the HO in terms 
of the triplet d-density state (addressed below) to explain the observed anisotropic 
magnetic susceptibility.E^HUl While this is an interesting idea, so far it has not been 
able to provide a mechanism for superconductivity, which must be related to the 
HO state. 

An early attempt, with some family resemblance to the skyrmionnc mecha¬ 
nism discussed here, is the non-BCS mechanism of superconductivity suggested by 
Wiegmann,!^ as an extension of Frohlich mechanism in one dimension to higher 
dimensions. The crux was the concept of spectral flow. Consider gapped fermions 
in a static potential. Assume that the chemical potential fjL lies in the gap. When we 
change the potential adiabatically, the energy levels of the fermions move around. 
Typically the levels cannot cross but there are potentials such that an adiabatic 
and smooth variation creates even number of unoccupied states below /i, or forces 
some occupied levels to cross this level. For this spectral flow to occur, the variation 
of the potential, a soliton, must necessarily be topological. The index theorem then 
relates the topological charge of the soliton to the number of levels crossed. This 
phenomenon produces a compressible liquid, a superconductor. 

More recently, several interesting papers have led to discussions of supercon¬ 
ductivity in single and bilayer graphenes. Grover and SenthiP^ have provided a 
mechanism in which electrons hopping on a honeycomb lattice can lead to a charge- 
2e skyrmionic condensate, possibly relevant to single layer graphene. As to bilayer 
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graphene, a charge-4e skyrmionic condensate has been suggested by Lu and Her- 
buPl and Moonfi^ See also the earlier work in Ref. of charge-e skyrmions in a 
quantum Hall ferromagnet. 

The starting point is layered condensed matter systems with weak interplanar 
tunneling. For realistic application to URu 2 Si 2 these planes are along the body di¬ 
agonals of a body-centered tetragonal crystal containing the active U atoms. This 
obscures the main features of the skyrmionic mechanism greatly. In addition, both 
crystal field effect and spin-orbit coupling must be taken into account, as was pre¬ 
sented in the previous papers,!^® although, in the end, they turn out not to be of 
primary importance for the mechanism of superconductivity discussed here. Also, 
the nesting vector appropriate for this material should be (0, 0 , 2tt/c), where c is the 
height of the unit cell.® In order to reveal the principal aspects, in the present paper 
we shall simplify by considering the weakly coupled layers to be in the XT-planes 
and the nesting vectors to be (tt/o, tt/o, 0) and its symmetry complements, where 
a is the spacing of the square planar lattice, which we shall set to unity. Only a 
minimal band structure involving nearest and next-nearest neighbor hopping will 
be incorporated. We shall eschew all other complicating details that can be found 
in our published papers.®® 

There are two points that are crucial to our work. The first is rather simple: in 
the density wave state considered here, there are also Goldstone modes that can 
be easily seen by integrating out the fermions resulting in a non-linear a-model 
involving a unit vector N, the form of which is entirely determined by symmetry 
and is given below. The Goldstones are spectators to charge-2e skyrmions that can 
possibly Bose condense, or their fractionalized avatars that lead to a paired BCS 
state at T = 0. At finite temperatures, however, they could lead to interesting 
behavior.l^ The second point is more subtle: we assume that the hedgehog configu¬ 
rations are absent. This would require that the energy of the skyrmions be smaller 
than the single electron gap, a question that is likely to be model dependent. If this 
assumption is correct, the transition from the st-DDW state (discussed in Section 
2) to the superconducting state will correspond to a deconfined quantum critical 
point ai T — 0, which otherwise would be a hrst order transition, as in Landau the- 
oryU-OlMim In other words, skyrmion number conservation in each plane is crucial 
to the fractionalized mechanism explored here. 

We beg the readers’ patience in reading this manuscript. The path through this 
review is as follows: In Section 2 we begin with the relatively ill-understood density 
wave states of higher angular angular momentum. A case of crucial importance 
is what we term as the st-DDW. The state is a topological spin Hall insulator 
with quantized spin Hall conductance. In Section 3 we show that the quantum 
fluctuations from the mean field state define skyrmionic texture of charge 2e, zero 
spin and zero angular momentum. In principle, this charge 2e skyrmions can Bose 
condense, akin to the seminal work of Wiegmanrfi^ (and even earlier by Frbhlich). 
Although this in itself is interesting, we push the argument further in Section 4 
and show that fractionalization of skyrmions can result in chiral superconductivity, 
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which appears to have been observed URu 2 Si 2 . Section 5 then discusses the endpoint 
of our tortuous path. A pesky notational issue is the vector notation: if no confusions 
arise we will set Q to be simply Q. Similarly, wherever possible we shall set h = 
c = e = 1. 


2. Density Wave States of Non-Zero Angular Momentum 


Density wave states in correlated electron systems can be defined by the angular 
momentum quantum number and the fundamental nature of the condensatesA 
superconductor is a condensate of Cooper pairs, that is, the condensation is in the 
particle-particle channel. Thus, the antisymmetry of the wave function provides a 
strict restriction on the spin function. If the orbital function is symmetric, the spin 
function must be antisymmetric and vice versa. In contrast, the density wave states 
are condensates of bound pairs of electrons and holes. Because there are no require¬ 
ments mandated by exchange between the two distinct particles, the orbital wave 
function cannot constrain the spin wave function. Although angular momentum is 
not a strict quantum number in a crystal, we will continue to use it as a metaphor— 
the proper classification is in terms of the symmetries of the point group. 

For a superconductor (. = 0,1,2,... dehne s-wave (spin singlet), p-wave (spin 
triplet), and d-wave (spin singlet) condensates, etc. For a particle-hole condensate, 
that is, a density wave, £ = 0 comes in two varieties: the spin singlet version dehnes 
the familiar charge density wave (CDW) and the triplet version the spin density wave 
(SDW). The £ = 1 comes also in two versions and involves bond order. The case 
£ = 2, spin singlet, is not a wave of density at all, but corresponds to a staggered 
pattern of circulating charge currents, dubbed the d-density wave (DDW). The 
£ = 2, spin triplet, corresponds to a staggered pattern of circulating spin currents. 
The two-fold commensurate DDW breaks translation, time reversal, parity, and 
a rotation by 7r/2, while the product of any two symmetries is preserved. More 
specihcally, the DDW order parameter, of period-2, is dehned by 


4-|-Q,QCfe,/3) = (cosfca: - cosky) 6 a, 0 , 


( 1 ) 


where Q = (7r,7r). Note the similarity of the form factor with the d-wave supercon¬ 
ductor (DSC), but the factor of i signihes the breaking of time reversal symmetry. 
The Kronecker Sa^p reflects the fact that the order parameter transforms as identity 
in the spin space, hence a singlet. As mentioned above, the order parameter in the 
real space corresponds to a staggered pattern of circulating charge currents shown 
in Fig.[^ It has been proposecP^HUl p)Dyy gap is proportional to pseudogap 

T* in a cuprate superconductor. A triplet DDW order on the other hand is given 

by 


t ^ 

(A+Q,a^k,p) = i-^{cOSkx - COS ky) N ■ da,p, 


( 2 ) 


which clearly transforms as a triplet under rotation in spin space; ct’s denote the 
standard Pauli matrices. In the real space this order parameter corresponds to 
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Fig. 1. Staggered pattern of charge currents reflecting DDW order in a square planar lattice. 


circulating staggered spin currents. The unit vector N defines the direction of the 
spin. 

It has been known that triplet icjdx^-y^ order parameter corresponds to stag¬ 
gered circulating spin currents around a square plaquette.^^ The oppositely aligned 
spins circulate in opposite directions, as shown in Fig. This reminds us of topo¬ 
logical insulators where oppositely aligned edge-spins travel in opposite directions. 
However, there is no topological protection because the bulk is not gapped, but 
is a semimetal instead. A more interesting case is the st-DDW order parameter 
{iadx^_y 2 + dxy), where a = ±1 for up and down spins, with the quantization axis 
along z. In the momentum space, the order parameter will be 

and the form factors are 

Wo 

Wk = ^(cos k^-cosky), (4) 

Afc = Aq sin/ca; sin fcy, (5) 

corresponding to the d,^ 2 _y 2 and d^y components, respectively. The current pattern 
is unchanged but the hopping matrix elements along the diagonals are modulated 
by Afc, as shown in Fig. 

The singlet chiral {idx^_y 2 + d^y) density wave, however, breaks macroscopic 
time reversal symmetry and was employed to deduce possible polar Kerr effect and 
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Fig. 2. Triplet icrd^ 2 _y 2 density wave in the absence of an external magnetic field. The current 
pattern of each spin species on an elementary plaquette is shown. The state is a semimetal. On 
the other hand iad^ 2 _y 2 + dxy can be fully gapped for a range of chemical potential. An example 
is shown in Fig.[^ see Ref. 



Fig. 3. st-DDW order on a square lattice. The solid (open) arrowheads indicate the current 
directions of up- (down-) spins due to the triplet i<jd^ 2 _y 2 order. The solid lines indicate the 
nearest neighbor hopping. The dashed and dotted lines indicate different signs of modulations of 
the next-nearest neighbor hopping due to the singlet dxy order. 


anomalous Nernst effecPH^ in the pseudogap phase of the cuprates. 


2.1. Mean field theory of quantized spin Hall effect in st-DDW 

Unlike the semimetallic DDW, st-DDW has a non-vanishing quantized spin Hall 
effect for a range of values of the chemical potential. This is in fact a topological 
Mott insulator!^ because it is the electron-electron interaction that is necessary for 
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it to be realized. 

In the absence of an external magnetic field, the Hamiltonian is 
'Hd±rd -rN = Y^ 'i’lAk'ilk, 


( 6 ) 


where the summation is over the reduced Brillouin Zone (RBZ) bounded by ky ± 
kx = ±7r, and the spinor, 'I'^, is defined as (c^ A J 


cl, ^ ci, I, c^_(_Q ^). The chemical 


potential is subtracted for convenience, N being the number of particles.The matrix 
Ak is 


7lr = 


/ ek- R Afc + iWk 0 0 

Afe — iWk ek+Q — R 0 0 

0 0 €k-iL Ak- iWk 

\ 0 0 Ak + iWk ik+Q — fl J 

with a generic set of band parameters. 


( 7 ) 


Cfc = eifc + £2* (8) 

eifc = —2t(cosA:a; + cosfcy), e 2 fc = 4t'cos cos (9) 

Here e 2 k is possible next nearest neighbor hopping. Each of the two 2x2 blocks 
can be written in terms of two component spinors, 'tpk.a = {ck,(j,Ck+Q,af", where 
tr = ±1 = (t, i); for example, for the up spin block we have the eigenvalues (± 
refers to the upper and the lower bands respectively) 

Afe,± = t2k — R A Ek, Ek = yje\k + + A^. (10) 

and are plotted in Fig. see Ref. |21 Since up and down spin components are 
decoupled, the Chern number for each component can be computed separately. Note 
that while (e 2 fc — /i) is present in the eigenvalues, it cannot enter the eigenvectors, 
because the identity operator commutes with the Pauli matrices. After diagonalizing 
the Hamiltonian, we can obtain the eigenvectors 

where (0(a:) is the step function) 


( 11 ) 


2 1 /I _l_ \ 

Ok = arctan(^) + 7r0(-Afe). 

Ak 

The Berry curvature, Ho-,± is 

^a,± = iVk X (^l.±(k)IVfcl®a.±(k)) 

From the eigenstates, the Berry curvature can be written as 

Ok^ 


Vt 


cr,± 


= iVk X [(m± - u±)Vfe{icr —)]. 


( 12 ) 

(13) 

(14) 

(15) 

(16) 
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Since u±, v±, and 9k only depend on kx and ky, only the z component, is 

non-zero, which is given by 




cr d 


±cr 


‘^El 


■ilk. 

d9k 

d 


dky 

dk. 

< 

Wk 

ilk 


dWk 

deik 

dkx 

Sfcx 

dkx 

dAk 

dWk 

deik 

dky 

dky 

dky 


eifc 

Ek 


dkx 


(17) 


From the above determinant, we can see that the Berry curvature will be zero if one 
of Afc and Wk is zero, so we need a mixing of dx 2 -y 2 and dxy to have a non-trivial 
topological invariant. If we define a unit vector ho- = ho/lhol, where 

ho = (Afc, -aWk, eifc), (18) 


the Berry curvature can be written as 

^(T,± = ^2^'^ 

More explicitly, the Chern numbers are 

d?k, 


Na,± = 


IRBZ 

= ±a 


27r 




( 7 ,± 


^ dkx ^ dky 


(19) 


[ (sin^ h + si^^ cos^ ky) 

Jrbz 2tt hj. 


( 20 ) 


= ±( 7 . 


We can focus on the lower band as long as there is a gap between the upper and 
the lower bands. Then, 


N = = 0 

Aspin = - A;,_ = (-1) - 1 = -2 


( 21 ) 

( 22 ) 


irrespective of the dimensionful parameters. Note, however, that the Chern numbers 
vanish unless both Aq and IFo are non-vanishing. The quantization holds for a range 
of chemical potential /i within the gap. 

For the fully gapped case, the ratio of the dimensions of the quantized spin Hall 
conductance to the quantized Hall conductance should be the same as the ratio of 
the spin to the charge carried by a particle, since in two dimensions both quantities 
have the scale dependence that is. 



So, the quantized spin Hall conductance will be 


(23) 


a 


spin _ 

xy ~ 


h 

h 2e 


N ■ = 

’ snin 


e 

27r 


(24) 
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Fig. 4. Energy spectra, Xk^±+fj,, corresponding to st-DDW. Here, for illustration, we have chosen 
Wo = t and Aq = —t and for illustration the band parameters are t = O.lSeV, t' = 0.3f. For the 
chemical potential, fi, anywhere within the spectral gap, the lower band is exactly a half-filled and 
the system is a Mott insulator, unlike the semimetallic DDW at half-filling. 


As long as time-reversal symmetry is preserved, we will still have Kramers degener¬ 
acy in our Hartree-Fock state, and therefore the edge modes resulting from topology 
will remain protected. We will dub this state as quantum spin-Hall insulator (QSHI). 


3. Low Energy Effective Action 

We now go beyond mean field theory. For simplicity consider the Hamiltonian (ad¬ 
dition of longer ranged hopping will not change our conclusion^ 


^ = E V'l 

k,a,/3 


<5o/3(T"eifc + r-Afe) - (a ■ N)^0TyWk 


'4’kp, 


(25) 


where the summation is over the RBZ. Here r® (i = x,y,z) are Pauli matrices 
acting on the two-component spinor. We expand around the points Ki = (f, f) 
and K 2 = (—§,§), the two distinct nodal points in the absence of the d^y term, 
and K 3 = ( 0 , 7 r), the nodal point in the absence of the dx 2 -y 2 term. This allows 
us to develop an effective low energy theory by separating the fast modes from the 
slow modes. 

After that we make a sequence of transformations for simplicity: (1) transform 
the Hamiltonian to the real space, which allows us to formulate the skyrmion prob¬ 
lem; (2) perform a 7 r /2 rotation along the r^'-direction, which allows us to match to 
the notation of Ref.l^ for the convenience of the reader; (3) label ipKi+q.a by V'zct) 
since Ki is now a redundant notation; (4) construct the imaginary time effective 
action, with the definition ip = Finally, after suppressing the spin indices, 

and with the definitions 7 ° = t^, 7 ^ = r^, and 7 *^ = —r^, we obtain the effective 
action in a more compact notation: 
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S = / cPx ipj 

1 = 1 , 2 “' 

.lEn 


- i^°dr - 2itj^{r]jd^ + dy) 


+ i —{a ■ N)j'^{-r]jd^ + dy) +ir]jA q 


(fxtp3 


- i-f°dr - Wo{a ■ N)'y^ 


tp3, 


where 771 = 1 and 772 = —1. There is no spatial derivative in the ips term since the 
expansion around the point K 3 = ( 0 , 7 r) is 


_ ^0 79 



(26) 


where the second (and higher) order derivative terms are irrelevant operators. The 
first term behaves as a mass term at the K 3 point. 


3.1. The charge and spin of a skyrmion 


First we will compute the charge of the skyrmions in the system.^i^ Consider the 
action around Ki = (^, ^) when the order parameter is uniform (say, N = z). The 
results for K 2 = (—f, f) and K 3 = ( 0 , 7 r) follow identically. We have shown above 
that in this case the non-trivial topology leads to a quantized spin Hall conductance 
in st-DDWstat^as long as the system is fully gapped. The spin quantum Hall effect 
implies that the external gauge fields and couple to charge and spin currents, 
respectively. Then, 


= / d^x'ipi 


- i^°dr + 7°(^r + Y^r) ” {d. + Oy) 


+ 2t^%Al + + AI+ ^-^AD + *^a^Y(-a, + dy) (27) 


Wo (7^ 

---a: +AI + -a;) + zAo 


i’l- 


The non-vanishing transverse spin conductance implies that the low energy effective 
action obtained from integrating out the fermions for the gauge fields is given by 

S3,eS=^Jd^xe^^^^A';,d^Al (28) 

and the charge current is induced by the spin gauge field 

f, = (29) 

Consider now a static configuration of the N field with unit Pontryagin index 
in the polar coordinate (r, 0 ): 


N{r, 9) = [sin a{r) cos 9, sin a{r) sin 9, cos a{r)\ 


(30) 
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with the boundary conditions a(r = 0) = 0 and a{r -A oo) = tt. This field configu¬ 
ration associates with a Skyrmion, and now the action is 

• VFo, 


Si = / 


- dr - f2t7 “idx + dy)+i — {d-N)-/^{-dj: + dy) + iAg 


V'l (31) 


We can perform a unitary transformation at all points in space such that 

C/1■(CT•iV)^7 = CT^ 

and dehne ip = and ■0 = 'ijj'W. Plugging into the above equation, we obtain 


.IPo 


- i"f^dr - 2itj^{dx + dy) +i — a^'y^{-dx + dy) + iAq 


S'! = y d^xijj'i 

+ [ d^xi;[ -ij°{U'^drU)-2it-f^{U'^d^U + U^dyU) 




+ i^a^l^{-U^d^U + U^dyU) 


To proceed, we write down the explicit form for U{r, 9), which is 

C a(r) • ct(r) —jO \ 

cos -Y -sm-^e \ 

• a(r-) i(r) 

sm^^e cos j 

and after applying the boundary conditions as r —i oo, we have 

17(r^oo,6l) = 

Therefore, in the far held limit, we have 

t/1'(r^oo,6l)(9,,l7(r^oo,d) = ^0 

isin0, _ 

= ( - 

r 


(32) 


(33) 


(34) 


U^r 


, 9 )dyU{r 


,0) = 


—e 


0 e \ .cos^ „ , / 0 —e 


0 


, i cos 9, ^ 

= (—K 


Putting into Eq.(|32[), now we have 

„o 


Si = J d^xtjj'i 

+ f d^x^'i 


.LFn 


(35) 


-Ydr - 2iY{dr + dy) -I- i^cf''l'^{-dx -f dy) + iAq 




2i7"(/. + fy) + - fy) 


V'i 


(36) 


where fy = —iWdyU. Equating the above equation and Eq. (27), as r —l oo, we 


obtain the explicit form of the gauge helds in the new basis in the far held limit: 

(37) 


4 c _ .c .s _ 2sin0 _ 2cos0 
- - Y ' y~ r ■ 
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In other words, the process of tuning the order parameter from to [a ■ N{r,9)\ 
in the original basis is equivalent in the new basis to adding an external spin gauge 
held 

- 2 sin 6* ^ 2 cos 0 „ 

A® =- X H- y 

r r 

= (38) 

r 

A Skyrmion with unit Pontryagin index in the + dxy state induces a 

spin gauge held A® = ‘^6. The total hux of this gauge held is 


$" = y d^xz • V X A" 

rdO • (-) = 47r. (39) 

r 

Suppose we adiabatically construct the Skyrmion conhguration N(r, 6) from the 
ground state z in a very large time period Tp oo. During the process, we effectively 
thread a spin gauge hux of 47r adiabatically into the system. The transverse spin Hall 
conductance implies that a radial current will be induced by the dir spin gauge 
hux of A®(t), which is now time-dependent: A®(t = 0) = 0 and A®(t = Tp) = A^. 
That is. 



m = -^dtAl{t). (40) 

As a result, charge will be transferred from the center to the boundary, and the total 
charge transferred during the process can be computed by performing an integral 
over the boundary and time: 




dt / rd0j^{t) 


= -j\de^[AUrp)-AU0)] 

= -^{2n) = -2. (41) 

Therefore, after putting back the unit of charge e to the expression, we obtain 
a Skyrmion with charge 2e. The result is identical for the nodal point K 2 . Since 
U^drU = 0, Sz terms do not couple to the gauge helds generated by the Skyrmionic 
texture. As a result, adding Sz terms does not affect the gauge helds obtained from 
Si and S '2 terms: 


A'^ = 0; A* = -e. (42) 

r 

Therefore, we have our hnal result: a Skyrmion in the st-DDW system carries charge 
2e. 

One can also verify the adiabatic result by a diherent method by computing the 
Chern numbers.l^ The charge and spin of the skyrmions are associated with the 
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coefficients of the Chern-Simons terms by the following relations: Qskym 
and S'skyrmion = <^11, where Ci and C 2 are 


Cl = 


C 2 = 


^fiu\ 

247r^ 

247r2 


Tr 


Tr 




dk„ 


dK 

-1 


dkx 

y-1 




dG 


-1 


dkx 


= Cae 


(43) 


(44) 


where G is the matrix Green’s function and the trace is taken over the spin index 
cr and other discrete indices. 

If the Green’s function matrix is diagonal in the spin index, then the Chern- 
Simons coefficients for up and down spins can be computed separately. 


V(G,) = [I d^G. 


9G-1 ac-i dG7^^ 


dkf^ dky dkx 


(45) 


(46) 


and Cl =N{G^) +N{Gi), C 2 = A/’(G^) — A/’(G^). It can be shown that 

G~^ = ioji — T - ha 

with ha being the Anderson’s pseudospin vectoil^ of the Hamiltonian, where for 
i(jdj.-i_y -2 -\-dxy system, we have ha dehned above. The Chern-Simons coefficient for 
spin cr can be written as 


d^k 


47r 




dh„ 


(47) 


where ha = ha/\ha\ is a unit vector . Here Ci and C 2 are the total Chern number 
N and the spin Chern number Agpin defined in the previous section, respectively. 
Explicitly, Ci = —1 -I- 1 = 0 and C 2 = —1 — 1 = —2; thus the results are the same 
as above. 

Because a Skyrmion in the system carries integer spin, it obeys bosonic statistics 
and may undergo Bose-Einstein condensate. As a result, the charge-2e Skyrmion 
condensate will lead to a superconducting phase transition. But what about its 
orbital angular momentum? It can be shown that it is zercP resulting in a s-wave 
singlet state. This is a bit surprising given the original d-wave form factor. 


3.2. The non-linear cr-model 

It is also possible to derive the corresponding non-linear cr-model, which also demon¬ 
strates the existence of skyrmions, as we have neglected the hedgehog contribution 
because skyrmionic excitations have lower energy than the particle-hole excitations 
required for tunneling between the layers. Finally we are ready to write down the 
non-linear cr—model. The effective action for the field N can be written as (correct¬ 
ing a mistake in Refl3]) 


5'eff « — 
91 
1 


dTd^x[\dxN\'^ + |9yiVp] + ~ y dTd^x\drN\ 


= - / dTd^x[\drN\^+vi {\d,N\^ + \dyN\^)] 


9 


(48) 

(49) 
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with the identification 


1 

9 


33 ^ 2(fc2 + W§y 

k 

V- ^0 

33 _ 2[fcg+£:^] 

3i 

2(fc^+WQ )2 


(50) 

(51) 


where = 4t^[A:a; + A;y]^ + ^[—kx + ky]'^ + Aq, and k = {ko,kx,ky). We can rescale 
the spatial coordinates by absorbing the coefficient and rewrite the effective 
action in a more compact form 


Ses « ^ J dT(fx\dfj,N\‘^ 

This completes our derivations of the NLcrM. 


(52) 


4. Chiral d-wave Pairing 

We have seen above that Bose condensation of charge 2e skyrmions could result in 
a charged superfluid with zero angular momentum and zero spin. For application to 
URu 2 Si 2 , it is necessary to find a mechanism for chiral d-wave superconductor EZEll 
This will be accomplished by fractionalization of skyrmions into merons and anti- 
merons. A meron is half a skyrmion. A mapping on the surface of a sphere of the 
two-dimensional plane covers only half the solid angle 27r instead of 47r. 

For the purpose of orientation, consider a phase diagram in which we introduce 
a quantum parameter A in addition to the parameters pressure, P, and temperature 
r, as shown in Fig. see. Ref. HI A controls Wo(A) such that Wo(A < Ac) = 0 and 
IFo(A > Ac) yf 0. Since an isolated meron costs logarithmically infinite energy for 
A > Ac, merons and antimerons appear as bound pairs in skyrmions. The length 
scale of the confinement potential grows when approaching the deconfined quan¬ 
tum critical point Ac, where it diverges. Therefore, the skyrmions fractionalize into 
merons and antimerons, because there is no confinement at that point. It is assumed 
that the hedgehog configurations are suppressed because the particle-hole excita¬ 
tions are of much higher energy.^ Therefore the skyrmion number is conserved in 
the two-dimensional AT-planes. 

The state at T = 0, P = 0 is connected, as is the entire superconducting state, by 
continuity. Ac is a deconfined quantum critical point. The suppression of hedgehog 
configurations is crucial to the existence of the deconfined quantum critical point. 
These fractional particles, merons and antimerons, are present at Ac, but not in 
either side of it. We presume that Ac can be computed from a suitable microscopic 
Hamiltonian; for instance, A may be a function of the on-site Coulomb interaction 
U, the nearest neighbor direct interaction E, and the exchange interaction J in an 
extended Hubbard model.l^ 
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Since merons have topological charge they should have zero overlap with band 
fermions and therefore cannot be expressed in terms of local band fermonic opera¬ 
tors. This is not unprecedented: the creation/annihilation of Laughlin quasiparticles 
in the fractional quantum Hall effect cannot be expressed as any local function of 
the band fermions. 


T 



A 


Fig. 5. The proposed phase diagranfll with the quantum parameter (A), pressure (P), and tem¬ 
perature (T) axes. Here A is a tuning parameter such that VFofA < Ac) = 0 and iyo(A > Ac) ^ 0. 
Ac is a deconfined quantum critical point between the QSHI and superconductivity as T = P = 0. 
Tjjo and Tc are the HO and superconducting transition temperatures as P = A = 0, respectively. 
Along the P axis, Pc indicates the phase transition between the HO and superconducting states, 
while Px indicates the phase transition between the HO and the large moment antiferromagnetic 
states Ellll In some literatures Pc coincides with Px^ which, however, does not affect our main con¬ 
clusion. In addition, there should be phase boundaries in the X-P and A-T planes, which are not 
the main purpose of this work. 


4.1. The interaction term 

Consider a skyrmion with a flux of Itt (in unit of hcle)^ as a composite of a meron 
with a flux of 27r and an anti-meron with a flux of —27r, as shown in Fig. We 
can imagine that at the critical point, the merons get deconfined. These fractional 
particles emerge as natural degrees of freedom right at the deconfined quantum 
critical point .1^ Thus, we can study the pairing instability which results from the 
interaction between these fractional particles. These deconfined merons only emerge 
Ac, but not on either side of it. In the st-DDW phase, the fractional particles are 
confined in skyrmions while in the superconducting phase they are bound into a 
Cooper pairs (Fig. see Ref. El). 

Let ^{r) be the creation operator of a meron at r, where s = ± labels the 
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Fig. 6. The merons ^{ 7 ^ and 'ip^_ '0+ creates a meron with iV(r ^ 0) = 2 and 

N{r —>■ 00 ) = ; ip'[_ ^(n) creates a meron with 7V(r ^ 0) = —2 and N{r ^ 00 ) = _ 

Each meron above is half the skyrmion: A composite of '0^ ^(f) and makes one skyrmion; 

see Ref. m 

flux of ±27r and cr =t(t) for up- (down-) spin. Here carries charge of ±e, so 

a skyrmion, a composite of A ^{f) and carries charge e — (—e) = 2e and 

flux 27r — (—27r) = Att. The subtraction is because destroying a particle with charge 
—e and flux — 27r is equivalent to creating a particle with charge +e and flux -f27r. 
Therefore, a pairing of {fpl a-i'^4’1' cr'(^) results in a charge 2e superconductivity 
for s = s'. For singlet pairing of the merons, a = —a'. 

Notice that the meron-anti-meron pair which constitute a skyrmion is not the 
same as the meron-meron pair which lead to the superconductivity. At the decon- 
fined quantum critical point, the merons are not bound in skyrmions, so they can 
interact with the merons within other skyrmions. This may be the reason why a 
skymion has zero angular momenturrP while a Cooper pair formed by the merons 
may have nonzero angular momentum. In other words, we are not studying the 
internal structure of a skyrmion. Instead, we are studying the pairing mechanism 
due to these fractional particles emerging at the deconflned quantum critical point. 

To formulate the interaction between the merons at the critical point, we imagine 
a situation where the gauge flux density V x A is produced by a meron. Suppose an¬ 
other meron passes the region of the created gauge flux, it will feel the Lorentz force. 
Because these particles carry charges, the motion of the merons also produces charge 
currents, and a magnetic interaction is induced between the two particles We 
also assume that the length scale of a meron is much smaller than the distance 
between merons, so that we can treat them as point particles. This is a critical 
assumption, as is the assumption that in the continuum limit discussed below, the 
merons or antimerons have kinetic energy /2m, with a finite mass. Going beyond 
these assumptions is outside the scope of this paper. 

The coupling between the gauge field and the charge current can be described 
by the imaginary time Lagrangian density, which is 


T^int — ^ 


(53) 
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Fig. 7. The deconfinement and pairing of the merons.^The up (down) arrows indicate a flux of 
27r(—27r). The solid (open) circles indicate merons (anti-merons). The colors are associated with 
the meron texture in Fig. [6| Left: in the st-DDW phase, a pair of a meron and an anti-meron is 
confined in a skyrmion. Middle: at the deconfined quantum critical point, the merons, which are 
fractions of a skyrmion, emerge. Right: In the superconducting phase, the merons are bound again 
into a Cooper pair. 


Therefore, the interaction Hamiltonian is 

= J d^rj{f) ■ A{f^ 

= (54) 

The gauge flux is induced by the spin texture, and adding a meron ijj± i® 
equivalent to threading a gauge flux of ±27r at f. 

So, the total gauge flux density and the density of the merons are related through 
the following relation, 

(v X d(f))^ = (55) 

S,(T 

where (V x A)z = dxAy — dyA^, and the charge density of the merons is 

(56) 

s,cr 

In addition, we have the continuity equation, which describes the relation between 
the density and the current of the merons, 

^ + V-JM=0 (57) 

With the continuity equation, and assuming that the kinetic energy of merons is 
e/c = ^ with the effective mass of merons m, we may write the current in terms of 
the V± <t(^ operator. 


(58) 
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The Fourier transform of Eq.(55) and Eq.(|58[) gives 




(59) 


and 




(60) 


Now the current is expressed in terms of '0^ ^ The next step is to write the gauge 
held in terms of '0^ ^ ^ as well. This can be done by choosing the Coulomb gauge, 
V • A[7^ = 0, or iq ■ Aq = 0. This means we can write Aq in the following form, 

2 

A = ^ A , (61) 

where q^ =ql + q^ and Aq = \Aq\ = ^A^^^+A^ y. 

Plugging into Eq.(59), we obtain the expression for Aq and hence Aq, 

^9 = ^ XI ^'^l,s,c^'>Pk+q,s,a {qyX - q,,y) . (62) 

k,S,(T 

From Eq.( |60[ ), we have the interaction Hamiltonian 

X X X k2 V’li+g,si,(Ti^fcl.Sl.fflV'fc2,S2.<T2^'=2+g,S2,<T2 


q,ki,k2 CTi,(T2 


q 


(63) 

Since we are interested in the Cooper pairing, we may set ki + k 2 + q = 0 and sum 
over q to obtain 

27ri ^ ^ —(fcl X fc2)z_,.t _,.t ./. 

^int — XXX *1^2 J \2 '^-k2,si,cn'^ki,si,ai'^Pk2,e2,<72'^-kl,S2,<y2 

kuk2^^S2cri,a2 I 1+ 2| 

Setting ki = k and ^2 = —k', si = S 2 = s, ai = —02 = cr, for singlet charge 2e 
pairing, the interaction Hamiltonian can be written as 

2 ikxk% i t 


iLmt — / , / , • 

k.k' 

Under time reversal, we have 


|fc- fc'l' 


'~'^k',s,ar-k>,s-cr^-k,s,-a^k,s,a 


(65) 


, 27ri (fc X P)^ I I 

= “~XX Tr- ^Wk',s,a'^-k\s,-a^-k,s,-a^k,s,a, 

- - \k — k'\^ 


k,k' 


^ Ui 


( 66 ) 


so the time reversal symmetry is broken in the interaction term, and we expect the 
superconducting order to break time reversal symmetry. 
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4.2. The BCS gap equation 

With the kinetic energy term, the total Hamiltonian is 

'Htotal = yk'k'4’\i ,s,i^-k,s,i'4’k,sA^ (6"^) 

fc,s,(7 k,k',s 


where — p, with the chemical potential p, and the summation over a was 

carried out. The BCS interaction is given by 


Airi (k X k')z 
m \k-k'\2 


( 68 ) 


Here we can see that the s = + and s = — parts are identical, so we will suppress 
the flux index s in what follows. In other words, we have 


^total — 'y ) ^k'4^j^ fj'4^k,cr T ^ ) ^k'k 






k'A^-k'A 




fe.T 


k,(T 


k,k' 


(69) 


With bk = {ilj-k,iipk,t) and bl = and defining, as usual, = 

— J2k' ^kk'bk', the superconducting gap, we get (note the distinction with the am¬ 
plitude of the dxy gap in st-DDW) 


by = 


(V'-fe'.tV’fc'.t) 

‘2Vey + |A^i 


i-2/(Ve + iA^fP) 


(70) 


where f{x) is the Fermi-Dirac distribution function. The BCS gap equation at T = 0 
is therefore 


A 


SC 

k 


Vkk' 


A 


sc 

hL 


2%/fr 


(71) 


where Vky = —Vyk- Below we shall also solve the gap equation at finite tempera¬ 
tures. 

Because of the imaginary prefactor of the potential, the s wave pairing would not 
be a solution since Al*^ cannot be real. The interaction is similar to the one discussed 
in the half-filled Landau level problerrP^ as well as the one in the context of the 
hole-doped cuprates.^^ However, there are some differences we would like to stress. 
In Ref. HOl the flux attached to a particle is ire instead of 27r to study fractional 
statistics by varying e. The particles in their system are spinless fermions, so they 
obtained a pairing state with odd-parity. The interaction is also different from the 
one in Ref. mi because we express p{'r) and A{'r) in terms of the ips,a{i^ operator 
differently. In Ref. HU the resulting interaction depends on the sign of s, so the 
s = ± part leads to a T id^y) superconductivity, respectively. As a result, 

the addition of these two components gives a dx 2 -y 2 superconductivity, but not a 
chiral state in cuprates. The merons and antimerons provide two identical copies of 
chiral superconductors. 
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4.3. The solution of the gap equation 

We begin with the ansatz for I wave pairing, 


(72) 


where 4 >k denotes the direction of the wave vector, and we will choose it to be the 
angle between k and k' for simplicity. Plugging the ansatz into the gap equation, 
we have 

p 2 tt 


iAri = - 


_ k^-\-k'‘ 


Anm 


k'dk' 


\^P\ 


d(j)- 


sin 


i/Cfe' + 7o Afcfc/ - cos 


(73) 


where Xkk' = 2 kk' ■ 

The angular integral can be computed by performing a contour integral in the 
complex plane. To do this, we set 2 = and get 




dz 


z 2Xz-z‘^-l’ 

where we have used dcj) = X = Xkk' for simplicity, and converted the (j) integral 

into a contour integral around the origin with unit radius 

(z2 - 


Ii{X) = <j) dz-^ 


\z-X-VX^){z-X + VX^y 

where the poles are at z = A ± VA^ — 1. Since z = X + VA^ — 1 > 1 for A > 1, 
it is not enclosed by the contour. Thus, only 2 ; = A — VA^ — 1 contributes to the 
integral, and we get 

Ii{X) = 2m{X - VA2 - 1)' (76) 

As A = Xkk' = , we have 

ask>k' 


^kk' \l ^\k' ^ “ 


(A) as k < k' 


So, 


Ii{Xkk') = 2.TTi X 


/ ,\i 

I — I 

\k ) 


a,s k > k' 
{^y ask<k' 
Putting back into the gap equation, we have 




fk' 
dk' I ^ 




VW+W¥\ 


dk' ( y, 


k'\Ap\ 


k'J VW+W\ 


(77) 


(78) 


(79) 


Here we can see that |A^°| (x k * as fc —)■ 00 and |A|'^| oc fc* as fc —>■ 0, so we take 
the ansatz 

i 


iAri = 


A17 ^ 


(^) ask< kp, 
i^y as k>kF, 


(80) 
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where A|f = | with the Fermi wave vector kp. Assuming k = kp and plugging 


Eq.(|80[) into Eq.(|79[), we obtain 

fkp 


\sc 

AF = ^ 
2m 


dk'- 


(&) 


21 


dk'- 


t'(t) 


-21 


+ yj (.Si - S,y+Air {li) 

( 81 ) 

In order to perform the numerical calculation, we need to make the variables di- 


-21 


mensionless. Let x = we have 

K.-c’ ' 


ASCU2 

^sc ^ 

2 m 


^ 21+1 


„1-2Z 


dx 




dx 


.21 




Dividing both sides by e_F and letting u = , we get 


^ 21+1 


A-21 


= dx 


0 i/ {x^ — I)^ + Ji 


dx- 


p2Z /, ^(^^2 _ l'j2 _j_ y2^-2l 

We can further simplify the equation by setting y = x^ and obtain 


1 

u = - 


dy 


uy 


2 Jo ij{y- 1)^ + v?y’^ 


^ 1 


dy 


uy 


Viy- 1)^ + u‘^y~’' 


Eq. (84) can be solved by iteration using Mathematica, 

nl „.Z 


_ I 

'^out — 77 


dy- 


UinU 


1 


dy- 


UinV 


-i 


(82) 

(83) 

(84) 

(85) 


270 ^/{y - 1)^ + uf^y^ ‘2-Ji i/(y - 1)2 + uf^y ' 

For the singlet spin state, only the even I pairing is allowed, so we have the 
singlet d wave paring. The superconducting gap will be 

Ar = |Ar|e2*^^ 

where 

0.406eF as fc > kp, 

0.406eF (it) as fc < kp, 

As k < kp, the gap is 

Al‘^ (X k^ (cos 2(()fe + i sin 20fe) 

= {kl - kl) + 2ikxky, ( 88 ) 

so it will be {d^^-y^ +idxy) superconducting state. 

At finite temperatures, the gap equation becomes 


iAri = 


( 86 ) 


(87) 


Af 




(89) 
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Similar to the zero temperature case, the solution for ^ wave pairing will be A^^{T) = 
|A^‘^(T)|e*^'^'', and the magnitude of the gap is now 


iAr(oi 


forfc<fc;^ 
Af (T) {Iff , for fc > kp 


(90) 


with temperature-dependent gap Ap{T) = |A|^(T)|. 

At a given temperature, the gap equation can be solved by iteration as above. 
The temperature dependence of the gap A|f(T) can be found for all non-zero 
here we show some dominant channels, I = 1,2,3,4, in Figure The zero temper¬ 
ature gap and superconducting transition temperature are listed in Table From 





Fig. 8. Tempertaure dependence of the pairing gap for the angular momentum channels i = 
1, 2,3, 4 from top to bottom. Only £ = 2,4 are relevant for singlet pairing. 


Table one can see that the ratio 2A|f(T = 0)/(fcBTc) ~ 4.390 for £ = 2, which 
is comparable to the experimental value 4.9 from the point contact spectroscopy.^^ 
Notice that both the theoretical and experimental values are larger than the BCS 
value ^ 3.52. 


Table 1. The zero temperature gaps and transition 
temperatures for the angular momentum channels, 
1 = 1 , 2 , 3 , 4 . 


i 

1 

2 

3 

4 

Aff (T = Q)/ep 

0.916 

0.406 

0.264 

0.197 


0.413 

0.185 

0.121 

0.090 

2Af{T = 0)/fcsTc 

4.436 

4.390 

4.364 

4.378 


5. Discussion 

The purpose of this review was to explore the role of skyrmions in condensed mat¬ 
ter systems and to point out how it can lead to a novel chiral superconducting 
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state via fractionalization into merons and antimerons. Although the framework 
was cast in the context of HO state in a heavy fermion model, it is the chiral-d wave 
superconductivity that was the emphasis. 

We have shown that the spin texture in the st-DDW leads to skyrmions, which 
can then fractionalize. This is not possible without the mixing of the singlet and 
the triplet components of the d-density wave. The skyrmions acquire charge 2e as 
a result of spectral flow induced by threading flux through the center of the sys¬ 
tem. Our work should be interesting because it has been pointed out that, based 
on the charge, the thermal transport and the specific heat measurements that the 
superconducting state in the URu 2 Si 2 may have the form oc kz{kx ± 
which is also a chiral d-wave superconductor. A more detailed comparison with ex¬ 
periments is beyond the scope of this review. As remarked in the Introduction, we 
have purposely avoided various complications to get to the essence of the skyrmion 
physics. However, we are greatly encouraged by the beautiful PKE measurements 
by Schemm et aZ.,1^ which directly probes breaking of breaking of TRS in the su¬ 
perconducting state. 

However, there are a number of unresolved issues: 


(1) In the PKE measurements, there is a pronounced anomaly at T* ^ 0.8 — IK 

within the superconducting phase.^ This may involve the £ = 4 channels below 
the temperature This explanation may be supported by the fact that 

the ratio of the superconducting transition temperatures for £ = 4 and 1 = 2 
channels, ~ 0.486, is comparable to the experimental ratio of the 

anomaly temperature to the superconducting transition temperature T*/Tc ~ 
0.533 — 0.667. It is tempting to suggest that the subdominant order in £ = 4 
channel is excited by the large laser frequency. This could be very similar to 
^He where the subdominant pairing in the /-wave channel is visible only in 
collective mode measurements EllEll It would be interesting to vary the laser 
frequency, if possible. 

(2) A less direct measurement of broken TRS in the superconducting state was 
recently presented in Ref. 45 On the other hand, a NMR experiment!^ finds 
conflicting results of broken TRS in the HO state itself. There appear to be 
no data below ^K. So, we do not know how the NMR signature of the broken 
TRS found in the HO state relates to that below the superconducting transition 
temperature below ^ l.SAT. observed in PKE. Clearly further NMR experiments 
will be helpful to settle this issue; see also Ref.[5]for PKE at higher temperatures 
and its interpretation. 

(3) One of the experimental signatures of the HO state is a specific heat jump 
AC/T « 270 mJ/mol-K^ at Thoi followed by an exponential behavior below 
ThO) which can be fitted with a gap of « II meV.I^ In Fig. the numerical 
calculation of the specific heat of the st-DDW state is presented, where the 
superconducting state is not included. Here an exponential behavior can be 
seen for T < Tjjoj consistent with experiments. The chosen parameters were 
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Wq = 14 meV and Aq = 13 meV. This choice of the parameters results in an 
overall st-DDW gap of 11 meV; see Eq. 10 The specific heat jump AC/T « 344 
mJ/mol-K^ at Tho- Note that we did not include the temperature dependence 
of the st-DDW gap, so the results in the region T ^ Tho are overestimated, 
but this does not change our principal conclusions. We also do not know how 
to properly model the low carrier concentration semimetallic behavior at low 
temperatures. It may well be that the linear specihc heat is due to two-level 
systems,^ which can be tested by improving the sample quality. There could 
also be alternate explanations in terms of thermal excitations of merons. A 
more pedestrian reason could be that the three-dimensional Fermi surface is 
only partially gapped. This can be seen if we rotate to the diagonal planes, as 
in our previous work.Sl We have purposely avoided this rotation in the present 
review so as not to obscure the primary mechanism. 


C , 

— (mJ/mol-X'^) 



Fig. 9. The calculated specific heat of the st-DDW state plotted as C/T vs T/Thq. The line 
above Tfjo is obtained from the experimental data for URu 2 Si 2 .l^ 


(4) We assumed the existence of deconfined merons and antimerons at a deconfined 
quantum critical point for which it was necessary that hedgehog configurations 
were suppressedEHl^l^ However, we do not know if such a deconfined quantum 
criticality actually happens, although it results in an exciting possibility. As 
emphasized above, a bound meron-antimeron pair which forms a skyrmion of 
charge 2e and its subsequent BEC condensation is clearly a different mechanism 
from the superconductivity from the BCS condensation of merons. If so, there 
could in principle be a BCS-BEC crossover or a transition between the two. It 
would be interesting to explore this avenue. 

(5) There may be other orders competing with the superconducting state. One 
possibility is a Skyrmion lattice, which may be understood by considering the 
repulsive interaction between the Skyrmions. It will be similar to the Abrikosov 
vortex lattice due to the vortex-vortex repulsion. If such an interaction exists, 
the skyrmion lattice can compete with the superconductivity. However, we do 
not know how to formulate this problem. Maybe a more general question is 
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whether there are any other competing or coexisting orders due to the non¬ 
trivial spin texture.. The crystallization of skyrmions has recently been con¬ 
firmed in neutron scattering studies of the three-dimensional helical magnets 
MnSP and Fei_^Co^Si.>^ The real space imaging of a skyrmion lattice has 
also been reported in the chiral magnet Feo.sCoo.sSi.l^ 

(6) We assume the length scale of a meron is much smaller than the distance be¬ 
tween two merons, so the merons can be treated as point particles; i.e. in the 
dilute meron approximation. What would happen if we remove this restriction? 
The stability of the superconducting state with respect to Coulomb repulsion 
between merons must be seriously addressed. The pairing instability should be 
examined including the Coulomb repulsion, 

He = ^ VqP-qPq, (91) 

q 

where Eg = — is the two dimensional Coulomb potential with the dielectric 
constant e^. However, we do not know how the stability against the meron- 
meron repulsion relates to the Skyrmion-Skyrmion repulsion which leads to the 
Skyrmion lattice. 

(7) In the continuum approximation, we assumed that the kinetic energy of a meron 

,2 

is Cfc = ^, but we did not justify it. Although, this is the simplest possible as¬ 
sumption, it needs justiheation, including the definition of the mass of a meron. 
Perhaps a more general question is whether the results will hold if we assume 
a different form for the kinetic energy of the merons. What is the extension to 
the lattice version? 

(8) When we constructed the effective Hamiltonian for merons, we set a chemical 
potential p for them. How is this p related to the chemical potential in the 
original electron system? Notice that in the original electron system the chemical 
potential needs to lie in the gap so that the adiabatic argument for charge- 
2e skyrmions holds. A speculative answer to this question is that the meron 
chemical potential is half the skyrmion gap, easily obtained from the non-linear 
tj-model. Alternately, the meron chemical potential can be inferred from the 
Fermi energy determined by comparing with experiments. 

(9) An extremely interesting question is the role of the Goldstones, as inferred from 
the non-linear cr-model. Are they merely spectators or do they have important 
physical consequences? The coupling of the gapless edge modes and the gapless 
Goldstone modes in the QSHI is an open question.l^ 
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